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Abstract 

We investigate the scaling behaviour of quenched Wilson twisted mass 
fermions at maximal twist applying two definitions of the critical mass. The 
first definition uses the vanishing of the pseudoscalar meson mass mps while 
the second employs the vanishing of the PC AC quark mass mpcAC- We 
confirm in both cases the expected 0{a) improvement. In addition, we show 
that the PCAC quark mass definition leads to substantially reduced 0(a 2 ) 
cut-off effects even when the pseudoscalar meson mass mps is as small as 270 
MeV. At a fixed value of mpg we perform continuum limits for the vector 
meson mass my and for the pseudoscalar decay constant /ps and discuss the 
renormalisation constant Zy of the vector current. 
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1 Introduction 



Knowing the scaling behaviour of a given lattice action is an important step towards 
performing a controlled continuum limit of physical observables. The important 
result that Wilson fermions with a twisted mass term added £Q lead, at full twist, 
to an automatic 0(a)-improvement of correlators j2] without the need of additional 
improvement terms makes it particularly interesting to test this formulation of lattice 
QCD with numerical simulations. In this paper we will use lattice spacings ranging 
from a = 0.05 fm to a = 0.17 fm and pseudoscalar meson masses from mpg = 1180 
MeV down to mps = 270 MeV in order to perform a detailed scaling test of Wilson 
twisted mass fermions. In this way, we extend a first scaling test [3] substantially. 
All our simulations presented here are done in the quenched approximation, see also 
|H HE H3 H] f° r further quenched results of Wilson twisted mass QCD. Let us mention 
that also full QCD simulations with this approach have already been performed and 
proved to be very useful in studying the phase structure of lattice QCD with Wilson 
fermions [HI El EH HH EE2! • 

A basic ingredient to have automatic 0(a)-improvement is the tuning of the bare 
quark mass to its critical value. In the language of twisted mass QCD (tmQCD) 
this procedure corresponds to fix the twist angle ton/2. In the present work, we will 
employ two definitions for the critical mass. The first is the point where the pseu- 
doscalar meson mass vanishes, the second, where the PCAC quark mass vanishes. 
In the following we will refer to the first situation as the "pion definition" and to 
the second situation as the "PCAC definition" of the critical point. Both definitions 
should lead to 0(a)-improvement, but they can induce very different 0(a 2 ) effects, 
in particular at small pseudoscalar meson masses. Indeed, in ref. |H EI] we reported 
that the pion definition can have substantial 0(a 2 ) effects which are amplified when 
the quark mass becomes small and violates the inequality fi > aA 2 (where \x is the 
bare quark mass, i.e. the parameter which provides at full twist a mass to the pseu- 
doscalar meson). On the other hand, when the PCAC definition of the critical mass 
is used, these particular kind of 0(a 2 ) cut-off effects are dramatically reduced as was 
demonstrated in 0[7|. The effect of reducing 0(a 2 ) artefacts has been theoretically 
studied in chiral perturbation theory in refs. [TH ITo] and put on a more general basis 
in ref. |16j . 

In this paper, we employ the pion and the PCAC definitions of the critical mass 
to study the scaling behaviour of the pseudoscalar meson decay constant and the 
vector meson mass. The renormalisation constant Zy of the vector current is also 
presented for three values of the lattice spacing. 

2 Wilson twisted mass fermions 

Wilson twisted mass fermions can be arranged to be 0(a) improved without em- 
ploying specific improvement terms The Wilson tmQCD action in the twisted 
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basis can be written as 

S[U,ip,ip} = a 4 ^2ip(x)(D w + m + ifi^v^ifjix) , (1) 



where the Wilson-Dirac operator Dw is given by 
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3 1 

= Yl 2 + " ° V ^ ] (2) 

and V M and V* denote the usual forward and backward derivatives. We refer to 
[3] for further unexplained notations. The definition of the critical mass m c will be 
discussed in detail in the next section. 

We extract pseudoscalar and vector meson masses from the correlation functions 
at full twist (m = m c ): 

C a P (x ) = a 3 J2(P a (x)P a m a = 1,2 (3) 



X 
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C a A (x ) = ^E^W^(°)> a = 1 ' 2 W 



k=l x 
k=X x 

where we consider the usual local bilinears P a = ip'js 1 ^-^' = ^Ikl^^ and 
The untwisted PCAC quark mass mpcAC can be extracted from the ratio 

m _ ExjMo^) pa W> „ - i a (a) 



Using the exact lattice PCVC relation 

(dffi(x)O(0)) = -2^ a \P\x)Om a = 1, 2 (7) 

(where <9* is the lattice backward derivative, V® is the point-splitted vector current 
and O is a local lattice operator) we can also compute the pseudoscalar meson 
decay constant at maximal twist without requiring any renormalisation constant 
(see [I3HE1B]): 

/ps = ^|<0|P a |PS>| a = 1,2, (8) 
mp s 

where mps is the charged pseudoscalar mass and \PS) denotes the corresponding 
pseudoscalar state. 

If one now is interested in using the local definition of the vector current = 
Vu has to be substituted with Z V V^ in eq. (J7J) and the relation is now valid 
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only up to order 0(a 2 ) lattice artifacts. This relation allows for a determination of 
Zy through 



3 Definition of the critical mass 

The Wilson tmQCD action of eq. $1} can be studied in the full parameter space 
(mo,/i). A special case arises, however, when m is tuned to its critical value m c . 
In such, and only in such a situation, all physical quantities are, or can easily be, 
0(a) improved. The critical mass, or alternatively the critical hopping parameter 
k c = (2am c + 8) _1 , has thus to be fixed in the actual simulation to achieve a twist 
angle of ir/2. 

One possible definition of k c is to perform at vanishing twisted mass parameter 
an extrapolation of (mps«) 2 — > 0. However, as we said in the introduction, it leads 
to the presence of large lattice artifacts which are amplified at small quark masses. 
A better definition of k c (which reduces this kind of 0(a 2 ) artifacts) can be obtained 
by first determining at fixed non-zero twisted mass parameter the value of K c (fia), 
where the PCAC quark mass of eq. (JHJ) vanishes. The non-zero value of the twisted 
mass allows a safe interpolation in this case. As a further step, an only short linear 
extrapolation of K c (fj,a) from small values of \xa to fia — yields a definition of k c 
which is expected |14 [ IT5 | IT6] to lead to 0(a 2 ) lattice artefacts which remain small 
even at small masses. In ref. [Sj we gave an example of such a computation of the 
critical mass and found that indeed the cut-off effects at small quark masses are 
substantially reduced when compared to results obtained from the pion definition 
of K c . 

Recently, a definition of maximal twist from parity conservation has also been 
investigated in |7j. There, the critical masses m c (/xa) were not extrapolated to 
fia = 0, but were used at the respective twisted mass parameter at which they were 
determined. 

4 Numerical results 

In this section we will provide a comparison of Wilson twisted mass results for the 
pseudoscalar meson mass and the pseudoscalar meson decay constant as obtained 
using the pion and the PCAC definitions of k c . Moreover we will present results 
for the vector meson mass and the renormalisation constant Zy as obtained using 
only the latter definition of k c . Our simulations are performed for 6 different values 
of the lattice spacing, 0.05 < a < 0.17 fm. The fermion matrix was inverted for 
a number of bare quark masses in a corresponding pseudoscalar meson mass range 




a = 1,2, 



(9) 



where <9 M is the symmetric lattice derivative. 
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13 


5.70 


5.85 


6.00 


6.10 


6.20 


6.45 


a (fm) 


0.171 


0.123 


0.093 


0.079 


0.068 


0.048 


r /a 


2.930 


4.067 


5.368 


6.324 


7.360 


10.458 


L/a 


12 


16 


16 


20 


24 


32 


T/a 


32 


32 


32 


40 


48 


64 




pion definition (ftJ? lon ) 


N 

* meas 


600 


378 


387 


300 


260 


182 


jd±a 


0.0070 


0.0050 


0.0038 


0.0032 


0.0028 


0.0020 




0.0139 


0.0100 


0.0076 


0.0064 


0.0055 


0.0039 


/i 3 a 


0.0278 


0.0200 


0.0151 


0.0128 


0.0111 




H4a 


0.0556 


0.0400 


0.0302 


0.0257 


0.0221 




fi 5 a 


0.0834 


0.0600 


0.0454 


0.0385 


0.0332 




/i 6 a 


0.1112 


0.0800 


0.0605 


0.0514 


0.0442 






0.1390 


0.1000 


0.0756 


0.0642 


0.0553 






PC AC definition (< CAC ) 


N 

meas 


600 


500 


400 




300 






0.0070 


0.0050 


0.0038 




0.0028 




H2CL 


0.0139 


0.0100 


0.0076 




0.0055 




/i 3 a 


0.0278 


0.0200 


0.0151 




0.0111 




\i^a 


0.0556 


0.0400 


0.0302 




0.0221 




fi 5 a 


0.0834 


0.0600 


0.0454 




0.0332 




fi^a 


0.1112 


0.0800 


0.0605 




0.0442 




fija 


0.1390 


0.1000 


0.0756 




0.0553 






0.0200 


0.0144 


0.0109 




0.0080 




flgd 


0.0420 


0.0302 


0.0228 




0.0166 





Table 1: Simulation parameters and statistics (N meas ) 

of 270 MeV < mpg < 1180 MeV using a multiple mass solver (MMS), which is 
explained in the appendix, on O(200) — O(600) (depending on the (3) gauge field 
configurations generated with the Wilson plaquette gauge action. See table [T] for 
an overview of the simulation points and table El for the values of n c obtained with 
both determinations. 

In table El we present the values of the pseudoscalar masses (in lattice units) for 
a set of 7 values of /ij (i = 1, . . . , 7) of the twisted mass, chosen at the different 
/3's in order to roughly give the same value of the pseudoscalar mass in physical 
units. For the PCAC definition of k c we added two more intermediate masses in the 
MMS, denoted as fis and fig in the tables. In table |U we present the corresponding 
values for /ps, as obtained from eq. (JSJ). The goal of this paper is to perform the 
continuum extrapolation of /pg at a fixed value of mps, for a number of values of 
mps- The first step of this procedure consists in interpolating the value of fps r o for 
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,,.pion 


K PCAC 


5.7 


0.169198(48) 


0.171013(160) 


5.85 


0.161662(17) 


0.162379(93) 


6.0 


0.156911(35) 


0.157409(72) 


6.1 


0.154876(10) 




6.2 


0.153199(16) 


0.153447(32) 


6.45 


0.150009(11) 





Table 2: Critical values of the hopping parameters obtained from the vanishing of the 
pseudoscalar meson mass (K% lon ) and from the vanishing of the PC AC mass (k^ cac ). 

the chosen values of mpsr . These values are close to the simulated ones, and so 
even a linear interpolation is usually sufficient. Since we do not want to extrapolate 
out of the range of simulated masses, the lowest mass that can be reached is given 
by the highest value of mpsr (in both sets of data) corresponding to /ii, which is 
given by the point at j3 = 6.45 with the pion definition of k c . This point corresponds 
to m PS = 297 MeV *. 

In fig ^ we plot fps r o as a function of (a/Vo) 2 for pseudoscalar masses in the 
range 297-1032 MeV, obtained with both definitions of k c . For values of beta large 
enough, the values /ps^o show, with both definitions of k c , a linear behaviour in 
(a/r ) 2 . This nicely demonstrates the 0(a) improvement for both definitions of 
k c . However, for the pion definition we notice that the effects in 0(a 2 ) are rather 
large, in particular at small pseudoscalar meson masses of 297 MeV and 377 MeV. 
In contrast, the PCAC definition reveals an almost flat behaviour as a function of 
(a/r ) 2 even at these small pseudoscalar meson masses. This confirms the results 
of j3] that the PCAC definition provides a better definition of the critical mass 
with substantially improved scaling properties of physical observables, especially at 
small quark masses. In order to take the continuum limit we identify the scaling 
region where the data are well described by corrections linear in (a/ro) 2 . This turns 
out to start at /3 — 6.0 with the pion definition of k c and at j3 = 5.85 with the 
PCAC definition of k c . The values in the continuum limit are obtained separately 
by performing linear fits to the data in these two regions. The results of these fits are 
shown in fig. ^together with the data. It is very reassuring that these independent 
linear fits lead to completely consistent continuum values. 

In view of the better scaling behaviour of the data obtained with the PCAC 
definition of k c , for which the continuum extrapolation is safely under control already 
with data corresponding to values of j3 in the range [5.85 — 6.2], we decided to use 
only these data to obtain the final results in the continuum limit. As already said, 
data with the PCAC definition of k c are affected by much smaller lattice artifacts at 
small masses. In particular the minimal pseudoscalar meson mass that corresponds 

*Throughout this work we use the value rp = 0.5 fm. 
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13 


5.70 


5.85 


6.00 


6.10 


6.20 


6.45 








mpsa 


(ft;? ion ) 

V C J 






0.2455(23) 


0.1682(26) 


0.1385(66) 


0.1129(41) 


0.1004(27) 


0.0720(28) 




0.3237(16) 


0.2256(22) 


0.1764(42) 


0.1482(27) 


0.1298(23) 


0.0914(27) 


fi 3 a 


0.4434(11) 


0.3122(19) 


0.2373(32) 


0.2030(21) 


0.1768(17) 




jjL^a 


0.6272(09) 


0.4452(14) 


0.3335(22) 


0.2865(15) 


0.2463(15) 




fi 5 a 


0.7767(09) 


0.5535(12) 


0.4134(17) 


0.3534(13) 


0.3037(13) 




H^a 


0.9074(09) 


0.6488(13) 


0.4839(16) 


0.4130(13) 


0.3546(12) 




\xia 


1.0255(08) 


0.7358(12) 


0.5491(14) 


0.4676(12) 


0.4021(11) 










mpqa I 

-TO \ 


K PCAC) 
C / 




filCL 


0.2323(18) 


0.1640(23) 


0.1217(66) 




0.0934(24) 




/J>2 a 


0.3245(15) 


0.2289(17) 


0.1708(50) 




0.1276(21) 




fi 3 a 


0.4598(12) 


0.3232(13) 


0.2396(33) 




0.1779(18) 




\i^a 


0.6564(10) 


0.4606(11) 


0.3403(22) 




0.2492(13) 




fi 5 a 


0.8114(10) 


0.5701(10) 


0.4214(17) 




0.3071(12) 




/i 6 a 


0.9451(10) 


0.6658(09) 


0.4925(14) 




0.3588(10) 




fija 


1.0647(09) 


0.7530(09) 


0.5579(14) 




0.4062(09) 




fi$a 


0.3892(14) 


0.2741(15) 


0.2038(40) 




0.1519(20) 




flgd 


0.5678(11) 


0.3984(12) 


0.2948(26) 




0.2160(16) 





Table 3: Pseudoscalar meson masses mpsa for all simulation points. 

to Hi is now 272 MeV. We thus choose nine pseudoscalar meson masses in the range 
272-1177 MeV and extrapolate the corresponding values of fps to the continuum 
limit. The results are presented in table El and fig. 121 The values of /ps show a 
linear behaviour down to pseudoscalar meson masses of about 270 MeV without 
signs of chiral logarithms which, for this quantity, should only appear in quenched 
chiral perturbation theory (q%PT) beyond one loop. In the same figure we also plot 
the continuum values as obtained from results of the ALPHA Collaboration with 
non-perturbatively improved Wilson fermions [T^]. Notice that for this action, due 
to the presence of exceptional configurations, the smallest pseudoscalar meson mass 
that could be simulated was above 550 MeV. In table El we also present the results of 
a linear extrapolation of our data to the chiral limit (performed on the six smallest 
masses) through which we compute the values of the pion and kaon decay constant 
f w and fx (the latter in the SU(3) symmetric limit). The ratio of the two gives 
fx/ ' fn — 1-11(4), which is 10% smaller than what is obtained experimentally. This 
is however consistent with what was observed in previous quenched calculations [20J. 

The vector meson mass shows, in comparison to /ps, larger statistical fluctuations 
at small masses. This fact, together with the considerably larger lattice artifacts 
which affect the values obtained with the pion definition of k c , makes the continuum 
extrapolation too difficult in this case. As a consequence we limit ourselves to 
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13 


5.70 


5.85 


6.00 


6.10 


6.20 


6.45 








J tr o \ 






ji\a 


0.0986(10) 


0.0782(13) 


0.0516(17) 


0.0466(14) 


0.0437(13) 


0.0329(13) 


/i2a 


0.1195(10) 


0.0890(12) 


0.0632(11) 


0.0546(09) 


0.0500(11) 


0.0361(11) 


/i 3 a 


0.1418(11) 


0.1003(12) 


0.0740(09) 


0.0623(08) 


0.0562(10) 




jjL^a 


0.1685(11) 


0.1149(12) 


0.0859(09) 


0.0716(08) 


0.0637(10) 




fi 5 a 


0.1902(11) 


0.1273(13) 


0.0949(09) 


0.0790(08) 


0.0698(09) 






0.2112(12) 


0.1390(14) 


0.1029(10) 


0.0858(08) 


0.0754(09) 




\xia 


0.2320(13) 


0.1501(14) 


0.1104(10) 


0.0919(09) 


0.0806(09) 










fpqa ( 


K PCAC) 
c / 




filCL 


0.1267(14) 


0.0894(14) 


0.0689(27) 




0.0512(16) 




/J>2 a 


0.1345(13) 


0.0947(13) 


0.0711(13) 




0.0532(13) 




fi 3 a 


0.1472(12) 


0.1025(12) 


0.0763(10) 




0.0567(10) 




\i^a 


0.1697(12) 


0.1159(12) 


0.0858(10) 




0.0633(08) 




fi 5 a 


0.1914(13) 


0.1284(11) 


0.0944(10) 




0.0694(08) 




/i 6 a 


0.2134(14) 


0.1402(11) 


0.1025(10) 




0.0751(08) 




fija 


0.2358(15) 


0.1518(11) 


0.1100(10) 




0.0803(08) 




fi$a 


0.1403(13) 


0.0983(12) 


0.0734(11) 




0.0548(11) 




flgd 


0.1589(12) 


0.1095(12) 


0.0813(10) 




0.0601(09) 





Table 4: Pseudoscalar meson decay constants fpsa for all simulation points. 

present only the data obtained with the PCAC definition of k c . The vector meson 
mass has been extracted from the correlators eqs. (J3|) and (0) with local source and 
Jacobi-smeared sink [21]. We observe that the tensor correlator systematically shows 
smaller statistical fluctuations and thus we report in table only results obtained 
from this correlator. 

In fig. El we show the results for the vector meson mass as function of (a/ro) 2 - 
Again, even for small pseudoscalar meson masses, the behaviour of the vector meson 
mass is almost flat in (a/r ) 2 , indicating that 0(a 2 ) lattice artefacts are also small 
for this quantity. We perform linear fits of these data as function of (a/ro) 2 which 
are represented by the lines in fig. El The continuum extrapolated values for the 
vector meson mass are presented in table |H1 and in fig. H As a function of the 
pseudoscalar meson mass squared they show a linear behaviour without signs of 
qxPT artefacts. In fig. 0] we also plot the continuum values obtained with non- 
perturbatively improved Wilson fermions [T9*] . 

In tableElthe results of a linear extrapolation (performed with the seven smallest 
masses) of our data to the chiral limit are given. Moreover, we use this extrapolation 
to compute the values of m p and of tuk* (the latter in the SU(3) symmetric limit). 
As already observed in quenched calculations where the scale is determined through 
To, these values turn out to be 10 — 15% larger than the experimental values. 
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13 


5.70 


5.85 


6.00 


6.20 






Vflya ( 


K PCAC ) 
c J 






0.716(67) 


0.589(43) 


0.458(30) 


0.306(26) 




0.773(36) 


0.591(19) 


0.451(20) 


0.317(20) 


/i 3 a 


0.854(20) 


0.628(09) 


0.467(12) 


0.339(11) 




0.973(15) 


0.701(05) 


0.517(07) 


0.378(05) 




1.076(11) 


0.765(04) 


0.560(06) 


0.415(03) 




1.178(09) 


0.834(03) 


0.614(04) 


0.452(02) 


Hid 


1.277(08) 


0.902(03) 


0.666(03) 


0.488(02) 




0.812(26) 


0.606(13) 


0.464(14) 


0.327(16) 




0.919(14) 


0.666(06) 


0.494(08) 


0.359(07) 



Table 5: Vector meson masses mya for all simulation points with the PCAC defini- 
tion of K c . 

As a last quantity, we computed the renormalisation constant Zy which is defined 
in eq. © and is to be obtained in the chiral limit. This renormalisation constant 
is expected to be equal, apart from lattice artifacts, to the one obtained for pure 
Wilson fermions. Also in this case we do not consider data obtained with the 
pion definition of k c . The reason is, as was shown in 0], that in this case the 
lattice artifacts affecting Zy become very large at small masses, thus preventing 
the possibility of performing a reliable extrapolation to the chiral limit. Using the 
PCAC definition of k c , lattice artifacts are well under control also at small quark 
masses (with the exception of (3 = 5.7 where we prefer not to perform the chiral 
extrapolation). In fig. El we provide one example of Zy as a function of the quark 
mass. It turns out that Zy is to a good approximation linear in (afi) 2 . In tabled 
we report the values of Zy extrapolated to the chiral limit together with the results 
from standard perturbation theory (SPT)[22] (where the parameter of the expansion 
is a^ PT = #o A^ 71 "); with g\ = 6/(3 the bare lattice coupling constant) and from 
boosted perturbation theory (BPT)j22| (where the parameter of the expansion is 
^bpt _ Q,SPTy-p^ with V = l/3(Re{TrP Atiy }) the average value of the plaquette). 
Non-perturbative determinations with Wilson fermions using different methods are 
present in the literature for (3 = 6.0 and give results in the range [0.57-0.74]. The 
spread is due to different 0(a) artifacts which affect the different determinations. 

5 Conclusions 

In this paper, we have explored the potential of Wilson twisted mass fermions to 
reach small quark masses and fine values of the lattice spacing. Using the PCAC 
definition of the critical mass the scaling region is found to start already at (3 = 5.85, 
for the observables investigated here and for masses down to 270 MeV. In fact, with 
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mps [GeV] 


fps [GeV] 


m v [GeV] 


0.272 


0.1500(66) 


0.904(102) 


0.372 


1538T45) 


937(53) 


0.432 


0.1572(40) 


0.955(42) 


0.514 


0.1631(36) 


0.978(33) 


0.624 


0.1724(33) 


1.027(26) 


0.728 


0.1823(31) 


1.083(19) 


0.900 


0.2002(29) 


1.198(12) 


1.051 


0.2161(28) 


1.313(09) 


1.177 


0.2283(28) 


1.410(07) 


0.0 


0.1439(37) 


0.876(25) 


0.137 


0.1452(37) 


0.883(25) 


0.495 


0.1617(45) 


0.973(27) 



Table 6: /pg and my in the continuum with the PC AC definition of k c . The values 
in the last three rows are obtained from a linear fit on the smallest 6 masses (7 in 
the case of my) and correspond to: the values in the chiral limit (first row); f n and 
m p (second row); fx and m^* in the SU(3) symmetric limit (third row). 





Zy 


yBPT 


yBPT 


5.85 


0.5982(4) 


0.69 


0.82 


6.0 


0.6424(4) 


0.71 


0.83 


6.2 


0.6814(3) 


0.73 


0.83 



Table 7: Zy from the PCVC relation (PC AC definition of K c ) and comparison with 
standard perturbation theory (SPT) and boosted perturbation theory (BPT). 

this definition of the critical mass, correlators, apart from being automatically 0(a)- 
improved, are only affected by 0(a 2 ) discretisation errors which remain small in the 
region of masses that satisfy the (order of magnitude) inequality fi > a 2 Is?Q CD . In 
the case of the pion definition the 0(a 2 ) effects are instead small only in the region 
H > ak 2 Q CD and, as observed in jl] and in the present work, may become quite 
relevant at masses in the range 250-450 MeV and for the range of lattice spacings 
simulated here. In this case the scaling region starts, for masses down to 270 MeV, 
around (3 = 6.0. However, at the two smallest quark masses, we included a point 
at (3 = 6.45 in order to be sure that we can safely control the continuum limit 
extrapolation. 

In the case of /pg we have explicitly checked that both definitions of the critical 
mass lead independently to consistent values in the continuum limit. For further 
simulations, the PCAC definition of n c is clearly preferable as it leads to considerably 
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smaller lattice artifacts at small quark masses, allowing for an enlargement of the 
scaling region. 

The results of this paper clearly reveal that Wilson twisted mass fermions allow 
for simulations at pseudoscalar meson masses of about 250 MeV without running 
into problems with exceptional small eigenvalues or large lattice artefacts. This 
statement holds, when the PCAC definition of the critical mass is used. This is a 
very important lesson also for dynamical simulations, a lesson that we could only 
learn through the detailed quenched study performed here. In the dynamical case, 
however, small pseudoscalar meson masses are very difficult to simulate and further 
development of algorithms is still crucial. For example, modern variants of the 
Hybrid Monte Carlo algorithm allow to simulate pseudoscalar meson masses of about 
380 MeV [21I2S1 without hitting the "Berlin Wall" [2H1I2Z1- We therefore believe that 
by using these algorithmic techniques together with Wilson twisted mass fermions 
at maximal twist, realised by tuning the twist angle through the PCAC definition 
of the critical mass, it becomes realistic to address the small quark mass region in 
full, dynamical lattice QCD. 
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Appendix 



Multiple mass solver for twisted mass fermions 

In this appendix we show that within the Wilson twisted mass fermion formulation 
it is possible to apply the multiple mass solver (MMS) [2H1 EH1 EOj method to the 
conjugate gradient (CG) algorithm. We will call this algorithm CG-M and give here 
the details of the implementation. 

The advantage of the MMS is that it allows the computation of the solution of 
the following linear system 

(A + a)x-b = (10) 

for several values of a simultaneously, using only as many matrix-vector operations 
as the solution of a single value of a requires. 

We want to invert the Wilson twisted mass operator at a certain twisted mass 
fio obtaining automatically all the solutions for other twisted masses Hk (with \nk\ > 
| /xo | ) - In the so called twisted basis the Wilson twisted mass operator is 

D tm = D + ifi kl5 r 3 , k = l,...,N m (11) 

where D is the standard massive Wilson operator D = D\y + mo (see eq. (0) and ((21)) 
and N m is the number of additional twisted masses. The operator can be splitted 

as 

Am = Am + ifak ~ ^0)75r 3 , Am = D + t/2 ol5 T 3 . (12) 

The trivial observation is that 

Am Am = DQdW + A ~ Jig , (13) 

where we have used 7s-Djy75 = D^ w . Now clearly we have a shifted linear system 
(A + a k )x — b = with A = Am Am an d &k = We describe now the CG-M 

algorithm in order to solve the problem {A + ak)x — b = 0. The lower index indicates 
the iteration steps of the solver, while the upper index k refers to the shifted problem 
with <7fc. 
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CG-M Algorithm 

4 = 0,r =p k = b,a-! = C-i = Co = 1,$ = Po = 
for % = 0,1,2, 



J 5 

frv, tv 
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a n /3 n (l - 4 s -) + a n _i(l - cr^ 

sn— 1 



fc Sn+1 
yi "n a/. 

X n+1 = X n~ >r a nPn 
%n+l %n OtnPn 

r n+ i =r n - a n Ap n 
convergence check 

n _ ( r n+l,r n +l) 

Pn+1 — / s 

vm "^nj 

p n +i = r n+1 + (3 n+ ip n 

/-k fyK 

ok a Sn+l^n 

Pn+1 — Pn+1 > fc 

Sn^n 

Pn+1 = Cn+l r n+l + Pn+lPn 

end for 



We give here the algorithm explicitly again, since it has a different definition of 
Cn+i compared to the one of [SHI- This version allows to avoid roundoff errors when 
a k — /4 ~ becomes too large. 

We remind that when using a MMS the eventual preconditioning has to retain 
the shifted structure of the linear system. This means for example that it is not 
compatible with even-odd preconditioning. 
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Figure 1: r /ps as a function of (a/r ) 2 using the pion definition (open symbols) and 
the PCAC definition (filled symbols) of the critical mass; fits are performed with a 
linear function in (a/Vo) 2 separately for each set. 
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Figure 2: Continuum behaviour of fps (filled circles) as a function of the pseu- 
doscalar mass squared (in physical units). We also plot the continuum limit obtained 
with non-perturbatively improved Wilson fermions (open squares) [T§] . 
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Figure 3: The vector meson mass as a function of (a/r ) 2 as obtained from the 
PCAC definition of the critical mass. We also show a linear fit in (a/r ) 2 of these 
quantity at fixed value of the charged pseudoscalar meson mass. 
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Figure 4: The continuum limit values of the vector meson mass (filled circles) as a 
function of the pseudoscalar meson mass squared. We also plot the continuum limit 
obtained with non-perturbatively improved Wilson fermions (open squares) |19j . 
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Figure 5: An example of the extrapolation of Zy to zero quark mass at (3 = 5.85. 
The data show a linear behaviour in (a[i) 2 . 
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